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ABSTRACT 

We detect variability of the Fe Ka emission line in the spectrum of X-ray pulsar 4U 
1626-67, correlated with changes in its luminosity and in the shape of its pulse proHle. 
Analysis of archival Chandra and RXTE observations revealed the presence of an 
intrinsically narrow Fe Kct emission line in the spectrum obtained during the source’s 
current high luminosity period. However, the line was not present during an XMM- 
Newton observation seven years earlier, when the source was ^three times fainter. The 
line is resolved by the high energy grating of Chandra at the 98% confidence level, 
and its small intrinsic width, a = 36.eV, suggests reflection off an accretion 
disk at the radius R « (T.SlgJ) x 10® cm assuming a Keplerian disk, viewed at an 
inclination angle of 20°. This value is consistent with the radius of the magnetosphere 
of the pulsar, suggesting that the line originates near the inner edge of a disk that 
is truncated by the magnetic field of the neutron star. Timing analysis of the XMM- 
Newton and RXTE data revealed a major change in the pulse profile of the source 
from a distinct double peaked shape during the high luminosity state when the line 
was present, to a much more complex multi-peak structure during the low luminosity 
state. We argue that the appearance of the line and the change in the shape of the pulse 
profile are correlated and are the result of a major change in the emission diagram of 
the accretion column, from a pencil-beam pattern at low luminosity, to a fan-beam 
pattern at high luminosity. 

Key words: Keywords from the MNRAS website 


1 INTRODUCTION 

4U 1626-67 is an accreting X-ray pulsar located at a distance 
of ~ 5 — 13kpc from the Sun llChakrabartvl [l9^ ). It has 
a pulsation per iod of Papin ~7.7 s and an orbital period o f 
Rorb ~42 min, (iMiddleditch et aDl98ll : [Chakrabartvlll998l) . 
Its short orbital period classifies it as an ultra compact X- 
ray binary (UCXB), a subgroup of low mass X-ray binaries 
(LMXBs) with orbital periods of less than one hour. The 
short orbital periods of UCXBs suggest such tight orbits 
that only an evolved compact donor can fit. UCXBs, most 
likely, consist of a white dwarf or a helium star that is ac¬ 
creting onto a neutron star due to Roche lobe overflow (e.g. 
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]h95; Verbunt fc van den Heuvell 19951: Delove fc Bildsten 
200.^ ; Delove. Bildsten. fc Nelemansl l200ij ). Due to their 


evolved nature, UCXB donors are expected to be hydrogen 
deficient. Depending on the initial conditions and the 
environment in which they are created (e.g. being part 
of a globular cluster) their donors can follow different 
evolutionary paths, leading to a variety of objects, rang¬ 
ing from non-degenerate He stars to C/O or O/Ne/Mg 


19861: Podsiadlowski. Raooaoort. & Pfahll 

2002; 

Yuneelson. Nelemans. & van den Heuveil 

2002; 

Bildsten & Delove| 2004|) . 


Since its dis covery dCiacconi et al.l Il972l : 
iRappaport et af] 1 19771 ) . 4U 1626-67 has been observed 
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by all major X-ray observatories. Furthermore, UV and 
optical spectra have been obtained by the Hubble Space 
Telescope and the Very Large Telesco pe, respectively. 
Analy sis of BeppoSAX observations by lOrlandini et al.l 
(Il998ll revealed the presence of cyclotron absorption lines, 
suggesting strong magnetic field of « 3 x 10^^ G, making 
the source the only strongly magnetised neutron star in 
an ultracompact binary known so far. Another intriguing 
feature of its X-ray spectrum is the presence of prominent 
Ne and O emission lines. First detected in ASCA data, 
analyzed by lAngelini et al.l (Il995l l , they have been studied 
extensively using high re solution spectroscopy (ISchulz et al.l 
l200ll : [Krauss et aI1l2007^ . The lines are believed to originate 
in C/O or O/Ne dominated hot plasma in the vicinity of 
the compact object, and suggest a donor whose chemical 
composition is enriched by products of later stages of 
nuclear burning, most likely, a C-O-Ne or O-Ne-Mg white 
dwarf. These hndings are further corroborated by the HST 
UV spectrum that revealed both emission and absorption 
features from C, O , and Si but lacked He emission lines 
dHomer et al.l [200^ 1 and by the VLT optical spectrum that 
featured prominent C and O emission li nes, but showed no 
evidence of H or He dWerner et al.llioo^ h 

The timing properties and the shape of the pulse pro- 
Hle of 4U 1626-67 have also been studied extensively since 
its discovery in 1977. Since then, 4U 1626-67 has tran¬ 
sitioned through steady episodes of both spin-down and 
spin-up phases and and is currently in a spin-up period. 
During its hrst spin-up period between 1977 and 1990 the 
source’s luminosity was estimated at « 10^^ erg/sec and 
its pul se prohle displayed a charac t eristic double pe aked 
shape dWhite. Swank fc Holtl Il983l : iBeri et ^ [20l3). As 
the source moved closer to its hrst observed torque re¬ 
versal, its pulse prohle shape started to gradually lose 
its double peaked shape (iBeri et al ][20ll. After the hrst 
torque reversal in 1990 t he source entered a spin- down pe¬ 
riod, its hux decreased (IChakrabartv et al.l 119971 1 and the 
pulse prohle changed to a broader shape tha t did not dis¬ 
play the previously observed distinct peaks (iKrauss et al.l 
I 2 OO 3 ). In 2008, 4U 1626-67 underwent a new torque rever¬ 
sal and e ntered a new spin-up pha s e that is ongoing un¬ 
til to day djain. Paul fc Duttal I 2 OI 0 I : ICamero-Arranz et al.l 
l2012fl . During the second spin-up period the source lumi¬ 
nosity increased by ~2-3 times and th e pulse prohle shape 
returned to the double peaked shape (I Jain^ P aul fc Duttal 
I 2 OI 0 I : ICamero-Arranz et al.l[201^ : iBeri et al.ll2014ll~ 

The pulse prohles of accreting X-ray pulsars exhibit a 
variety of shapes, ranging from simple sinusoidal-like pro¬ 
hles, to clear double-peaked shapes, and to more com¬ 
plicated broadened prohles with multiple peaks. The ob¬ 
served variety of pulses among diherent - or sometimes 
same sources, rehects a variety of possible emission pat¬ 
terns of the polar region of an accreting neutron star. 
Depending on the mass accretion rate (~ source lumi¬ 
nosity), the emission diagram of the accretion column 
may switch from a pencil-be am to a fan-beam pattern 
dBasko fc Sunvae^ Il975l . Il976lj . At low luminosities, be¬ 
low « 10^^ erg/sec, the high anisotropy of the photon- 
electron scattering cross-sections in a high magnetic held of 


the order of 10^^ G (ICanuto. Lodenauai fc RudermarJll97ll : 
iLodenauai et al] 11974 1. leads to formation of the pencil 
beam p attern of radiation, orien ted parallel to the accretion 
column dBasko fc Sunva^ll975l l. However, at high mass ac¬ 
cretion rates, corresponding to lu minosities above the crit¬ 
ical value of Lc ^ 10®^ erg/sec dBasko fc Sunvae^ Il976l : 
IWang fc Frankll98lll , a radiation dominated shock is formed 
at the distance of a ~few km above the neutron star surface. 
At high luminosities, the accretion funnel is hlled with high 
density plasma slowly sinking in the gravitational held of the 
neutron star, resulting in increased opacity in the direction 
along the magnetic held axis. Consequently, the emerging 
X-ray photons predominantly escape from the - optically 
thin - sides of the accretion funnel and the fan beam pat¬ 
tern of radiation is formed (Fig. [3)). Pencil-beamed emission 
is usually associated with single-pulse prohles or more com¬ 
plex shapes if one includes gravita tional effe cts, and diher¬ 
ent obscuration mechanisms (e.g. iMeszaro^ 1992). Double 
peaked prohles are indica tive of an emission pattern that is 
fan-beam dominated ( e .g. Nage il98ll:IWhite Swa.nk fc H^ 
1 19831: iPaul et al.lllg^. ll997l : iRea et al.ll2004ll . 

In the present paper we investigate the X-ray spectrum 
and pulse prohle of 4U 1626-67 during two diherent lumi¬ 
nosity states. We present results of spectroscopic and tim¬ 
ing analysis of an XMM-Newton observation performed in 
2003 “ during the spin-down low luminosity period, and the 
latest simultaneous Chandra and RXTE observations per¬ 
formed in 2010, during the current high luminosity, spin-up 
period. In Section [2] we describe details of our data extrac¬ 
tion along with their spectral and timing analysis, followed 
by interpretation of our results and discussion in Section [3] 
and conclusions in Section |4] 


2 OBSERVATIONS, DATA ANALYSIS AND 
RESULTS 

4U 1626-67 has been observed multiple times by all ma¬ 
jor X-ray telescopes. For the present work we focus on the 
XMM-Newton observation performed in August 2003, the 
latest Chandra HETGs observation performed in January 
2010 and the RXTE observation that was performed si¬ 
multaneously with the Chandra observation. The details of 
the observations used in our analysis are listed in table [T] 
The temporal resolution of the EPIC pn instrument aboard 
XMM-Newton was sufficient for timing analysis of the 2003 
data. For the 2010 data, simultaneous RXTE observation 
allowed for high quality timing analysis of the source’s light 
curve during the Chandra grating observation. Spectral anal¬ 
ysis was carried ou t using the XS PEC spectral fitting pack¬ 
age, version 12.8.2 (lArnaudlll996l ~l. Timing analysis was per¬ 
formed using the standard tools of the XRONOS timing anal¬ 
ysis software package, version 5.22. 


2.1 XMM-Newton 2003 observation 

During the XMM-Newton observation MOSl detector was 
operating in timing mode, while MOS2 and pn detec¬ 
tors were operating in imaging mode. The MOS2 detector 
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Table 1. Details of observations of 4U 1626—67 analyzed in this 
paper 


Instrument 

obsID 

Date 

Duration^ (ks) 

XMM-Newton 

0152620101 

2003-08-20 

58 

Chandra 

11058 

2010-01-14 

77 

RXTE 

P95338-05-01-00 

2010-01-14 

10 


^Duration of filtered observations. 


showed evidence of pile up. For this reason and since the 
effective area of pn at « 7keV is approximately five times 
higher, we only use the pn data for the present analysis. 


2.1.1 Spectral extraction and analysis 


We extracted the source spectrum from a 30" circle cen¬ 
tered at the source. Background was extracted in compli¬ 
ance with the latest EPIC calibration note^ from a source- 
free region at the same RAWY position as the source re¬ 
gion. Spectral extraction was performed using the standard 
tools provided by the XMM-Newton Data Analysis software 
SAS, version 13.5.0. The resulting spectrum was re-binned 
to ensure a minimum of 25 counts per energy channel. 4U 
1626-67 is known for displaying strong emission features in 
the 0.5-1.5keV range. The study of the low energy part of 
the spectrum of this source is beyond t he scope of this pa - 
per and has already bee n performed bv ISchulz et akl (1200 il l 
and iKrauss et ah 1(H). Therefore, in order to simplify our 
analysis, we ignored energy channels below 1.5 keV. The re¬ 
maining channels are sufficient to constrain the source con¬ 
tinuum around the iron line. 

The spectral continuum was fit with an absorbed black 
body plus power law model. The temperature of the black 
body was 0.30±0.01 keV and the power law photon index 
0.70 ±0.01. Our fit also required an exponential cutoff at 
7.911 q 23 keV with a folding energy of 29.2lg4keV. The 
spectral shape at this energy resembles a mild break in the 
power law, rather than a genuine exponential cutoff. Ad¬ 
dition of the exponential cutoff model improves our fit by 
a Ax^ of 30 for 2 dof. Best fit parameters are presented 
in table [2] and a re in agreement with the findings of previ- 


,_ ' ( f g- 

iKrauss et al.l l2007tl for the spin-down era of 4U 1626-67. 
We do not detect the iron line, placing an upper limit of 
2.4 eV at 90% confidence, for the equivalent width (EW) of 
a Gaussian emission line centered at 6.4 keV with a 36 eV 
width. The choice of the values for the line centroid and 
width was motivated by the corresponding best fit values 
in the 2010 data, which are described in section The 
1.5-10 keV luminosity, calculated from the XMM-Newton 
fit, is ~ 1.3 X 10^ ^ erg/sec, assuming a distance of 9kpc 
llChakrabartvll99^ . Luminosity in the 0.01-100 keV, extrap¬ 
olated from the best fit model, is « 8.4 x 10^® erg/sec. The 
data-to-model ratio vs energy is presented in Fig. [T] 


ous authors (e.g. Angelini et al.l[l995l : lOrlandini et al.l| 19981 : 



li _I_I_ I _I_I 

° 3 5 8 


Energy (keV) 

Figure 1. Ratio of the data to the continuum model for the 2003 
XMM and 2010 Chandra and RXTE observations. The data have 
been rebinned for clarity; the 3-8 keV energy range is shown. 


2.1.2 Light curve extraction and analysis 

For our timing analysis we extracted a 2-12 keV light curve 
from the pn data using standard SAS tools. Photon arrival 
times were corrected to the solar system barycentre. Using 
the tools provided in the XRONOS package version 5.22, 
we determined the pulse period during the XMM-Newton 
observation and created the pulse profile of the source in 
the 2-12 keV range. The pulse period is measured at « 
7.67547(2) sec and the corresponding pulse profile is pre¬ 
sented in Fig. O 


2.2 Chandra and RXTE 2010 observation 

The January 2010 Chandra observation of 4U 1627-67 was 
performed with the high energy grating. During the Chan¬ 
dra observation a simultaneous 10 ks observation was also 
performed by RXTE. For our spectral analysis we used the 
data obtained by both Chandra and RXTE. The RXTE data 
were also used for timing analysis. 


2.2.1 Spectral extraction and analysis 

We extracted the spectra of both the medium energy grat¬ 
ing (MEG) and the high energy grating (HEG). Extraction 
was executed using the standard toolqj provided by the lat¬ 
est CIAO software (vers. 4.6.1). In order to remain consis¬ 
tent with the XMM-Newton data analysis, we chose to ig¬ 
nore energy channels bellow 1.5 keV in the Chandra data as 


^ http://xmm2.esac.esa.int/docs/documents/CAL-TN-0018.pdf 


^ http://cxc.harvard.edu/ciao/threads/pointlike/ 
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well. We also ignored MEG energy channels above 5keV. 
The Chandra data were not regrouped and were fitted us¬ 
ing the stan dard method and t he weighting technique 
suggested bv IChurazov et ^ dlQQStl . We also extracted the 
source spectrum from RXTE-PCA standard-2 data. Spec¬ 
tral extraction and background subtraction were performed 
using the standard routines provided by the FTOOLS pack¬ 
age, following the guidelines described in the RXTE Guest 
Observer Facilitjj^. 

Our analysis of the 2010 Chandra spectrum shows 
that since the 2003 XMM-Newton observation the pho¬ 
ton index of the power law component has increased to 
a value of 1.13±0.03 and the black body temperature to 
0.43±0.01 keV. Most importantly, our analysis of the Chan¬ 
dra spectrum reveals the presence of an iron Kq line that 
was not present in the XMM-Newton observation of 2003. 
The line - modeled using a simple Gaussian - is located at 
S.SOIq qj keV, has a width of cr = 36.eV and an EW 
of 18l5 g Addition of the Gaussian line improves our fit 
by a 5x^ of 13 for 3 dof, giving the detection of the feature 
a 2.8 cr significance. 

The presence of the line is confirmed by our fit of the 
3-20 keV, RXTE-PCA data, which strongly required a line 
with the best fit centroid energy of O.Gdlg’j® keV and the 
width of cr = 40l^g° eV. The EW of the line in the RXTE 
fit, has a value of 35.ll^g geV and is consistent -- within 
1(7 errors - with the Chandra value. While RXTE lacks 
the spectral resolution to constrain the line with the ac¬ 
curacy of Chandra, its larger effective area yields a spec¬ 
trum with substantially higher signal-to-noise ratio, thus 
providing a detection of the line with more than dcr sig¬ 
nificance ((5x^ = 38.1 for 3 dof). The data-to-model ratios 
for the Chandra and RXTE observations are presented in 
Fig. [T] The plots are for the best-fit model of the contin¬ 
uum, without the Gaussian emission line. The presence of 
the iron emission line can be seen clearly, particularly in 
the RXTE data that have higher signal-to-noise ratio. The 
1.5-10 keV luminosity, calculated from the Chandra fit, is 
« 4 X 10^® erg/sec. Luminosity in the 0.01-100 keV, extrap¬ 
olated from the best fit model, is « 3 x 10®^ erg/sec. Best fit 
values for the spectral parameters obtained from Chandra 
and RXTE spectral analysis are presented in tabled 


2.2.2 Light curve extraction and analysis 

For the timing analysis of the 2010 observation we made use 
of the RXTE-PCA good xenon data that have a time reso¬ 
lution of 1 ps. Using the standard tools described in the pre¬ 
vious sections, we extracted and analyzed the source’s light 
curve. We measured the pulsation period at « 7.67797(6) sec 
and produced the corresponding pulse profile in the 2-12 keV 
range. The pulse profile is presented in Fig. [21 


® http://heasarc.gsfc.nasa.gov/docs/xte/recipes/pca_spectra.html 





Phase 


Figure 2. The 2-12 keV pulse profile of 4U 1626-67 from XMM 
and RXTE observations. The switch to a distinct double horned 
shape in the 2010 RXTE observation is evident. The profiles have 
been arbitrarily shifted in phase so that the minimum appears at 
the pulse phase 1.0. 

3 DISCUSSION 

We have performed spectral and timing analysis of two dif¬ 
ferent observations of 4U 1626-67 taken during low and high 
luminosity periods of the source. We detected a faint, narrow 
Kq emission line of iron at 6.4 keV in the 2010 Chandra and 
RXTE data, during the high luminosity state of the source. 
This feature was not present in 2003 (Fig.[T]) when the source 
luminosity was less than 10®^ erg/sec. Furthermore, the ap¬ 
pearance of the line coincides with a major alteration in the 
source’s pulse profile. During the 2010 observation - in which 
the iron line is detected - the pulse profile of the source has 
a characteristic double peaked shape that is radically differ¬ 
ent from the pulse profile of the 2003 observation. Below, we 
discuss the details of the iron emission line and propose that 
its appearance is caused by the modification of the emission 
diagram of the accretion column caused by the change of the 
mass accretion rate. The same modification of the emission 
diagram leads to changes in the shape of the pulse profile. 

3.1 The iron Kq line 

Parameters of the iron lines observed in the spectra of ac¬ 
creting X-ray pulsars can vary significantly from source to 
source. Depending on the environment and the ionization 
state of iron, the energy of the K-shell line can vary from 
« 6.4 keV in the case of fluorescent line of neutral and weakly 
ionized iron to « 6.9 keV in the case of the resonant line 
of highly ionized, hydrogen-like FeXXVI. The observed line 
width also varies in the broad range from a few tens of eV to 
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Table 2. Best fit parameters for the XMM-Newton, Chandra and RXTE/PCA 
spectra. The errors are la. 


Model parameter 

XMM-Newton 2003 

Chandra 2010 

RXTE-PCA 2010 

nH“ (10®® cm-®) 

1.00 

1.00 

1.00 

Power Law 

T 

0.70*’ ±0.01 

1.13T0.03 

1.04 ±0.02 

norm‘s 

5.81T0.07 

34.9T1.5 

14 7+0-5 

Black Body 

kT (keV) 

0.30T0.01 

0.43T0.01 

0.61±0.07 

norm^ 

1.19T0.04 

6.95T0.41 

^•O0_o.38 

Iron Line 

Centroid E (keV) 

6.39'= 

6.39l°;°? 

6 fi 4 + 0-20 

D-04_q 25 

Width a{eY) 

36.4'’ 

36.41®®;® 

<400 

FluxI 

<0.39 

7.66l°-i® 

7 cq+3.73 

i •«JO_3 g4 

EW (eV) 

<2.41 

18.o1®;2 

Qtc 1+17.4 

OO.i_ig 9 

-^1.5-10 keV 

1.29 

4.00 

2 . 02 *“ 

X®/dof 

1731/1693 

3073/3008 

38.15/38 

“ Parameter frozen at total galactic HI column density provided by the HEASARC 

nH tool dDickev Kr. Lockmanlll990l:iKalberla et alJl20n5r). 



^ With a high energy cut off at 7.91 ^q' 23 keV with an e-folding energy of 
29.2^g’3g keV, modeled using XSPEC model highecut. 

10~® phkeV“^ cm“^ s“^. 

L 39 /DJ 3 j^p^ were L 39 is luminosity in units of 10®® erg/sec and Dioicp^ is 
distance in units of lOkpc. 

Parameter frozen. 
f 10“® phcm~®s“l. 

® 10 ®® erg/sec. 

^ Calculated in the 3-20 keV range. 


> 0.5 keV. Typical values of the equi valent width of the line 
are of the order of 0.1 to 0.4 keV (e.g. White. Swank fc Hold 
1 198,31 : lOottwald et al.l 1 199,51 : iTorreidn et al.ll20iri l. From a 
theoretical point of view, the iron line in X-ray pulsars 
may appear due to reflection of the primary X-ray emis- 
sion off the accretion disk and/ or the donor star (e.g. 
iBasko. Sunvaev fc Titarchutd [l973 l. It can also b e due to 
Alfve n shell emission, as proposed bv iBas ^ lll980ll . In their 
work, iBas argued that the iron emission line originating 
in the moderately optically thick, highly ionised plasma of 
the Alfven shell are expected to be centered at > 6.5 keV, 
are broader and brighter, with equivalent widths exceeding 

> 0.3 keV. On the contrary, fluorescent iron line produced in 
highly magnetised X-ray pulsars by reflection from the ac¬ 
cretion disk and the donor star surface, are expected to be 
generally narrower, centered at «6.4keV and have a mod¬ 
erate equivalent width of the order of 50-100 eV. 


The iron Ka line detected in the spectrum of 4U 1626-67 
during the 2010 Chandra observation is centered at 6.4 keV, 
has a small width of «36eV and an the equivalent width 
of «18eV. These line properties are more consistent with 
reflection from either the cool surface of the donor star or 
from the accretion disk. Distinguishing between the two lat¬ 
ter possibilities, we note the following. For a binary system 


consi sting of a l.dMpi N S accretor and a 0.02Mq WD donor 
(e.g. IChakrabartvl [l9^ 1 . the Roche lobe of the WD sub¬ 
tends a solid angle of ~ 5.1 x 10“^ sr as viewed from the 
neutron star. Assuming that the emission source is located 
in the accretion column, at ~15 km above the the disk plane 
and that the disk is truncated at the magnetospheric radius 
(see eq. [T] and discussion below), the disk will subtend a 
solid angle of ~ 1.3 x 10~^ sr as viewed from the emission 
source, ignoring possible flaring of the disk. This is smaller 
than the solid angle subtended by the WD donor, however, 
disk flaring and non-isotropic emission pattern can change 
this number and the relative contributions of the disk and 
the surface of the donor star. Indeed, the width of the line 
measured by the high energy grating, « 36 eV, corresponds 
to velocities of « 1700 km/s. Such velocities are more typi¬ 
cal for the accretion disk, rather than for the surface of the 
donor star, suggesting that the majority of the line emission 
originates from the disk. Furthermore, actual contributions 
of the accretion disk and the donor star to the observed fluo¬ 
rescent line flux are determined by not only the solid angle, 
but also depend on the emission diagram of the primary 
emission and the angles of the line of sight to the normal of 
the surfaces of the disk and the star. The latter are mainly 
defined by the inclination of the binary system. As for the 
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former, in the case of the fan beam, for example, it is quite 
possible that much larger fraction of the accretion colnmn 
emission is intercepted by the accretion disk, rather than by 
the donor star. 

The magnetospheric radius of the neutron star can be 
estimated as follows fe.g. lOhosh. Pethick fc LamblllflT^ b 


\MV2GMnsJ 


( 1 ) 


The magnetic field of the neutron star in 4U1626-67 is 
know n from the cyclotron l ine measurements, B « 3 x 10^^ 
(e.g. lOrlandini et al.l [l998f l. Assuming a Mns = 1.4Mq 
neutron star of the radius of Rns = 12 km, accreting at 
the mass accretion rate of M = 2.7 x 10~® Mq/ yr corre¬ 
sponding to the observed luminosity of « 3 x 10^^ erg/sec, 
the magnetospheric radius is Rm ~ 6.8 x 10® cm. This 
value is very close to the pulsar corotational radius Rco = 
(G Mivs Bspin/4 7r^)^^® « 6.5 x 10® cm. Hence, the pulsar 
is expected to be near spin equilibrium (Bm ~ Rco), (e.g . 
IShakuralll975l : iLipunov fc Shakui^ll976l : IziolkowskilllOSSli . 
The accretion disk will be disrupted by the magnetic field 
close to the magnetospheric radius, i.e. to the first approx¬ 
imation, the inner disk radius can be es t imated as Rjn ~ 
6.8 X 10® cm (e.g. IPringle fc Reed 1 19721 : iRomanova et al.l 
l2012l 'l. The outer disk radius will e xtend to a substantial 
fraction (up to 0.8, e.g. iLaso^ 1200 il l of the Roche lobe ra¬ 
dius, which is ~ 2.2 x 10^® cm, assum ing a white dwarf donor 
of 0.02 Mq fe.g. IChakrabart^^ll998^f . 

The line is resolved by the Chandra HETG at the 
98% confidence interval. The 90% error interval for the line 
width, 25.1 — 51.7eV, corresponds to the range of line of 
sight velocities of « 1178 — 2426 km/sec. For the disk in¬ 
clination ang le of 20° (e.g. IVerbunt. Wiiers fc Burrnlll990l : 
IChakrabartvl [l9^ f . this corresponds to the 3D velocity 
range of 3443—7092 km/sec. Keplerian velocity reaches these 
values at the distance of (3.7— 15.7) x 10® cm from the com¬ 
pact object, with the best fit value of the line centroid - 
36.4 keV, corresponding to the distance of 7.5 x 10® cm. 
These numbers are consistent with our estimate of the mag¬ 
netospheric radius (Rm ~ 6.8 x 10® cm, see above), suggest¬ 
ing that the fluorescent line is produced near the magneto¬ 
spheric boundary of the neutron star, where the Keplerian 
accretion disk is expected to be truncated by the neutron 
star magnetic field. At these radii, an accretion disk that 
is he ated through viscus dissipation llShakura fc SunvaevI 
[Tizi), cannot reach a temperature higher than « 5 — 6x 
lO'^ K. Obviously, at these temperatures any iron that is 
present in the disk will be in a low ionization or neutral 
state and its fluorescent K- alpha line will be centered at 
«6.4keV fe.g. lBeardenlll96^ 1 as observed. 

Although the observed line energy and width, are con¬ 
sistent with what should be expected for this system, the line 
itself is quite faint. With an EW of «18eV, it is significantly 
fainter than the EWs of « 40-100 eV of iron lines usually ob- 
served in disk reflec tion spectra of non-pulsar LMXBs (e.g. 
ICackett et al.ll^OlOl 'l. This could be due to the fact that the 
truncated accretion disk around the highly magnetised neu¬ 
tron star subtends a significantly smaller solid angle than 
in the case of ’’normal” LMXBs. However, the non-isotropic 


emission diagram of the accretion column can compensate 
for the smaller solid angle, for example in the case of the fan- 
beam as discussed above (Fig. [S]). Indeed, narrow fluorescent 
lines of iron, of appreciable equivalent width in the ~ 70 eV 
range, were previously detected in the spectra of some ac¬ 
creting X-ray pulsars in low-mass X-ray binary systems. For 
example, a na rrow line at ~6.4 keV was detected in GX 1+4 
( EW^70eV. IPaul et al.ll2005ll and in Her X-1 (E W^65 eV 


lEndo. Nagase fc Miharall2000l : IXaik fc Paulll200al . The 6.4 
keV line, we found in 4U 1626-67, is significantly fainter than 
in those sources. 

If the reflected component is not strongly reduced 
due to the geometry of the s ystem (the detectio n of 
O and Ne emissio n features by iKrauss et al.l (l2007l ~l and 
ISchulz et al.l (1200il l suggests it is not), then the faintness 
of the iron line can be explained by the f act that 4U 
1626-67 is an UGXB wi th a G/O-rich donor (lSchuIz et_alJ 
200l|_ HomeretalJ 20021 : |NelemanSj_^(onkei_&_Steegh^ 2006: 


Werner et al.l 200611. In Kolio panos. Gilfanov fc BildstenI 


( 2013ll and Koliopanos et al.l ( 2014 1 we demonstrated that 
the iron Ka line is strongly attenuated in reflection spectra 
of G/O or O/Ne dominated disks. This is due to screening of 
the presence of iron by the overabundant oxygen. In partic¬ 
ular, in the G/O dominated material the main interaction 
process for an E ^ 7.1 keV photon is absorption by oxy¬ 
gen rather than by iron, contrary to the case of “standard” 
LMXBs with main sequence or red giant donors where the 
accretion disk has a chemical composition close to Solar. 

To furth er investigate this hypothesis, we ran ou r sim¬ 
ulation from iKoliopanos. Gilfanov fc BildsteiJ (l2013l l for a 
primary radiation with a power law spectrum with a pho¬ 
ton index of 1.1. For simplicity, we assumed that the primary 
radiation is emitted isotropically by a point source above 
the disk surface, in a lamppost configuration. We have col¬ 
lected the reflected emission for viewing angles in the range 
15° — 25°. We found that in order to produce an output 
spectrum with the iron Ka line with an EW of « 18 eV - as 
observed - we need an 0/Fe ratio in the disk that is 68 times 
the solar value. Note that the maximum value of the 0/Fe 
ratio that corresponds to the chemical composition of a C/0 
white dwarf - in which all hydrogen and helium has been 
converted to carbon and oxygen - is « 77 times the solar 
value. This result suggests that the donor star in 4U1626-67 
is a C/0 or O/Ne/Mg white dwarf, in a perfect agreement 
with the previous work which proposed that this source is 
a UGXB wit h a C/ O white dwarf donor (ISchulz et al.ll200ll : 
iKrauss et al.l [2007 1. To conclude, we note that a lamppost 
geometry is obviously not an accurate representation of the 
emission of the accretion column. However, it can serve as a 
sufficient first approximation which demonstrates that the 
small equivalent width of the iron line in the spectrum of 4U 
1626-67 can be easily explained in terms of reflection from 
a C/O rich disk. 


3.2 The origin of the iron line variability. 

The iron line was not detected in the 2003 XMM-Newton 
data, with the upper limit of « 2.4 eV (In) on its equiva¬ 
lent width. A possible explanation for the line variability is 
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Figure 3. Schematic representation of the pencil and fan-beam 
emission in an accreting X-ray pulsar. Note that the drawing 
is only aimed to illustrate the difference of the disk illumina¬ 
tion at different emission patterns of the accretion column and is 
not meant to realistically reproduce the geometry of an accret¬ 
ing strongly magnetised neutron star. In particular, for viewing 
clarity, the distance of the shock from the NS surface has been 
exaggerated, as well as the size of the neutron star itself, whereas 
the inner disk radius has been depreciated. 

suggested by the timing analysis of the source. In particular, 
we notice that the line is present during the high luminosity 
state of the source in which the pulse profile has switched 
to a characteristic double peaked shape. 

We argue that both the appearance of the fluorescent 
line of iron and the change of the shape of the pulse profile 
are the result of a major modification in the emission dia¬ 
gram of the accretion column. Such a modification may be 
caused b y the increase of t h e ma ss accretion rate, as pro¬ 
posed by iBasko &: SunvaevI (Il975l i. Observationally, in the 
case of 4U1626-67, the increase of the mass accretion rate 
manifests itself as a more than 3-fold increase of the source 
luminosity and acceleration of the neutron star spin. 

At low luminosity, the emission of the accretion column 
is concentrated in a beam that is o riented along the magneti c 
field axis (pencil-beam emission) iBasko fc SunvaevI (Il975ll . 
Depending on the angle between the rotation axis and mag¬ 
netic dipole of the neutron star, the magnetic field axis - 
and hence the beamed radiation - will be mostly directed 
away from the accretion disk (Fig. [3]). Due to the decreased 
flux towards the disk, disk reflection features in the observed 
spectrum will be significantly suppressed, including the iron 
line emission. Note that the latter is weakened further due to 
the C/O dominated chemical composition of the accretion 
disk. 

As mass accretion rate increases, the emission diagram 
of the accretion column changes to the fan-beam pattern 
IBasko fc Sunva^ (Il976l i. This modification in the emission 
diagram of the accretion column, causes most of the radia¬ 
tion to be beamed towards the accretion disk (Fig. [3]), re¬ 
sulting in a significant boost of the emission that is reflected 
off the disk and the appearance of the detectable Ka line of 
iron. As the disk is truncated by the magnetic field of the 


neutron star at large distance form the latter, relativistic 
broadening of the line is insignificant and the line is rather 
narrow. The line is also rather faint, due to the fact that 
the disk is made of C/O-enriched material. The fan-beam 
diagram of the emission from the accretion column further 
manifests itself throughout the characteristic double peaked 
shape of the pulse profile. 

In this picture, we can predict that the iron line should 
have been suppressed during the entire spin-down period 
from ~ 1990 to ~ 2008. In order to check this prediction, 
we analysed two earlier Chandra observations performed in 
2000 (with duration of 39 ksec) and 2003 (95 ksec), during 
the spin-down era. We detect no emission lines in the 6- 
7keV energy range in either of the these observations. We 
obtained an upper limit (1 cr) on the equivalent width of the 
6.4 keV iron line of 24 eV in the 2000 data and a 7.1eV in 
the 2003 data. In both cases the Chandra spectra were fitted 
with the same spectral model as XMM-Newton data, and 
produced similar values of best fit parameters. The 6.4 keV 
iron line was not detected in the 1994 ASCA data neither, 
although with a fairly unconstraining upper limit of of 33 eV 
(lAngelini et al.l 1199511 . In all three observations the source 
was in the same low luminosity regime. These non-detections 
of the line provide further support to the picture outline 
above. 


4 SUMMARY AND CONCLUSIONS 

We have discovered the appearance of a narrow Fe Ka emis¬ 
sion line in the spectrum of 4U 1626-67 correlated with a 
major change in the shape of its pulse profile. The line was 
detected at a moderate level of EW = ISl^ g eV during a 
period of high luminosity in 2010, when the pulse profile of 
the source had a characteristic double peaked shape. The 
line was not detected in 2003 with an upper limit of 2.4 
eV (Icr) when the source was in a low luminosity state and 
the shape of its pulse profile had a clearly different, broader 
shape. We argue that both changes are caused by the same 
reason -- a significant modification of the emission diagram 
of the accretion column, from a pencil beam to a fan beam 
pattern. This change was caused by an increase in the mass 
accretion r ate, a s was theoretically predicted by Basko & 
Sunyaev in Il976l . The lower than typical equivalent width 
of the iron line in the high luminosity state is due to C/O 
dominated chemical composition of the accretion disk. 
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